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An explicit relation between the dimensional loss (AD), entropy production and transport is 
established under thermal gradients, relating the microscopic and macroscopic behaviors of the 
system. The extensivity of AD in systems with bulk behavior follows from the relation. The 
maximum Lyapunov exponents in thermal equilibrium and AD in non-equilibrium depend on the 
choice of heat-baths, while their product is unique and macroscopic. Finite size corrections are also 
- - - ^ computed and all results are verified numerically. 
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Fractal structures in phase space have become the focus of attention in the understanding of the relationship between 
^ microscopic dynamics and macroscopic non-equilibrium physics [Q. In the escape-rate formalism, the properties of 
a fractal repellor are known to govern the transport ||^. In contrast, in boundary driven non-equilibrium steady 
. states, the stationary distribution is generally fractal, but the precise connection to transport is not fully understood 
' 101 . This reduction in dimension, AD, has been argued to be related to transport [^-||, although the only precise 
T— I understanding arises in the weak field limit of the Lorentz gas [|[, which does not have clear generalization to non- 
equilibrium many-body systems. The presence of fractals has also been used to demonstrate how the second law of 

Q' ', thermodynamics is consistent with time-reversal invariant, deterministic dynamics Q. Many of these issues require 
the understanding of the Lyapunov spectrum, whose analytic properties are known only in certain special cases [Q. 
In this letter we study general Hamiltonians coupled to two heat baths at different temperatures at opposite sides 
i of the system, generating heat flow. We will see that the type of heat baths chosen affect AD and the Lyapunov 
• ^ ' exponents, for the same boundary temperatures. We investigate the meaning and the origins of extensivity of AD 
(~| , in systems with bulk behavior, including finite size corrections, and relate them to transport. The new relations are 
verified numerically. We systematically study the system not only close to but also far from equilibrium, as well as 
^ the dependence on the heat-baths themselves. While the extensivity of dimensional loss has been disputed due to the 
^ incompatibility with local equilibrium we will see that this is not an issue. 

\ In systems with bulk behavior, to say that AD behaves "extensively" means that AD should remain relatively the 
T-H • same under the same local non-equilibrium conditions when we change the size of the system. For systems in thermal 
^ gradients, one might expect 

. ^ ' where the constant C" might depend on T . Such a behavior is particularly natural if one envisions that a continuum 
j-j limit of the theory ultimately exists. AD has been studied previously for color conduction sheared fluids 

pO|,^p|] and thermal conduction [Q , numerically. Analytic computations of the Lyapunov spectrum and AD have 
understandably been restricted to small or idealized systems such as the Lorentz gas . The physical properties 

are far from trivial; even whether AD generally arises has been an issue ||l2j. Extensivity of AD under thermal 
^ ' gradients has been analyzed jj], but the relation to transport and entropy production was not elucidated previously. 
5^ , Extensivity has been investigated in sheared fluids |^,^ and for color conductivity |^,^ . Study of the dependence on 
the number and types of thermostats or systematic analysis far from equilibrium have not been performed before. 

We first discuss the properties of AD, its extensivity and its relation to transport. Consider a system with cross- 
sectional area A placed in contact with two heat-baths at both ends having temperatures {Ti,T2). The Lyapunov 
spectrum ^^-^ 2 d distills the microscopic properties of the classical system and are the time averaged eigenvalues 
of the stability matrix (dipi/difj), where {^j}^^^ 2 d degrees of the freedom of the system, including 

the heat-baths. From the spectrum, we may compute the fractal dimension using the Kaplan- Yorke estimate Dky, 
which has the property that AD = D — Dky > in non-equilibrium and AD ~ in equilibrium systems. While 
the extensivity relation (|l|) is natural, the relation is ill-defined for two reasons: First, VT is in general not constant 
within the system and which VT we choose visibly affects the results |l3j. Second, it is unclear what D in the 
denominator should mean. It can be the number of degrees of freedom of the whole system including the thermostats, 
that without, or that only of the interior. These difficulties arise since Dky is a global quantity defined only for the 
whole system. An expression consistent with (|l|) in the near equilibrium limit that suffers no such ambiguity even far 
from equilibrium is (denoting the heat flux by J), 
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since J is constant throughout the system. Vi„ is the interior volume of the system. We shall derive this relation near 
equilibrium and further verify it numerically. 

Computing AI? requires the full Lyapunov spectrum, which involves evolving + D degrees of freedom. This 
apparently makes it difficult to explicitly check extensivity. Close to equilibrium, this difficulty can be overcome as 
follows Define ADmax.ADmin as 

AD = -Sk^ AD ■ = ^^^^ rsi 

where Xjnaxi\nin ^rc the maximum and minimum exponents. When AD < I, ADj^in = AD holds exactly. (More 
generally, when K — 1 < AD < K, we need only to compute X^j^i ^^'^ ^ lowest Lyapunov exponents to 

obtain AD.) This holds when the system is close to equilibrium. VT/T ~ f/L and AD - L"^'^ (in d-dimensions) 
so that AD is always small for large systems in 1-d. Since Xmin = —^max in equilibrium, ADmax should also be a 
good approximation to AD, close to equilibrium. We analyze these conditions systematically in the following and 
make these conditions more precise. ^jj the total rate of phase space contraction, can be computed from the 

equations of motion of {^Pj} alone and Xmax can be computed from evolving one tangent vector, so we need only to 
evolve 2D degrees of freedom to compute AD max — a huge reduction for large D. 

We now systematically investigate how the extensivity of AD arises. First note that the rate of entropy 

production, 5, for the system 0,^. We use this to derive a thermodynamic relation 

where k is the thermal conductivity. We have included the effects of boundary temperature jumps that behave as 
aJ when the jumps are not too big |p^ . a, which arises as a finite size correction, measures the efficacy of the 
heat-baths, which can be stochastic or deterministic, and can have significant effects as will be shown. Notice that 
(4) always holds, both close to and far from equilibrium and is independent of the type and number of thermostats 
used. Further, it is Vin that arises in the relation, rather than D. 

Xmax behaves as Xmax — -^maa; + ^{J^)^ where A^^^ is its equilibrium value. A^^^ is independent of Vm for large 
enough systems, but can depend on T and also on the thermostats used. Close to equilibrium, the above behavior 
of Y^f^i Xj and AJ^^^, when combined, explain the extensivity of AD as in (||). In this limit, the extensivity of AD 

arises from X^jLi and the thermostat dependence from Af^^^. Since Cd can be derived in the J limit, we 
derive 

1 2aK \ Vin — > oo 1 
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which relates macroscopic transport and entropy production to the microscopic AD. A subtlety needs to be mentioned: 
We have found that A^^j^. is consistent with having a finite large volume limit but have not proven this statement 
analytically. In fact, this difhcult issue is open and there have been conflicting results on the existence of the 
thermodynamic limit of X'^^^ |p^ . From our argument, we see that the extensive nature of AD requires that Af^^^, 
have a thermodynamic limit or vice versa. We note that in systems without a bulk limit, such as the FPU model, AD 
needs not be extensive and A^^^ needs not have a thermodynamic limit, which might explain some of the discrepancy 
seen in the previous literature |15t] . Far from equilibrium, AD and AD^aa; will be quite different, yet the extensivity 
relation (||) can and seems to still hold for AD. This reflects the deeper geometric significance of AD not present in 

ADmax • 

While the above theory is valid in any dimension, we now apply it to the 1-d cj)'^ theory described by the following 
Hamiltonian: 
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We choose the (j>'^ theory because it is a classic statistical model that naturally appears in broad physical contexts. 
Also, the statistical properties of the theory have been studied previously, including thermal transport which has 
bulk behavior [|l3|,^. We model the heat-baths dynamically by applying Nose-Hoover (NH) thermostats or their 
variants (demons) at the boundaries, as explained in |ll3|] . The interior includes the dynamics only of the (j)'^ theory. 
Temperature is defined unambiguously using the ideal gas thermometer, T{x) = (tt^). The number of (heat-bath) 
boundary sites thermostatted on each end, Nb, will be varied. We employ one set of thermostats per thermostatted 
site. We run the simulations with time steps of 10"'^ to 0.05 for 10^ to 10® steps, understanding its convergence 
properties and also checking that the results do not change with the step size. To obtain the Lyapunov spectrum, we 
use the method in . 

Let us first look at the dependence of AD, ADmax and ADmin with respect to J (Fig. Q). We see that close 
enough to equilibrium, all three quantities agree and display behavior as in (|) and (|). Remarkably, even far from 
equilibrium and well into the non-linear regime, AD has a robust behavior, but not ADmax,min- We further note 
that the boundary temperature jumps have a significant effect, the effect being larger for smaller T . For instance, 
using K = 2.83(4)r"i-35(2),a = 2.6(1) [||, 2aK = 300,6 for T = 0.1,2, respectively. 
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FIG. 1. AD (□), ADr„ax (o) and AD„ii„ (A) against - J in L = 11 (left), L = 41 (right) at T = 0.5 for the Nb = I case, 
behavior with and without finite size corrections are also shown. AD displays behavior even for large — J. 



Performing an analysis similar to those performed for Fig. |^, we can extract the proportionality constant AD/J^ 
for aparticular L, Nb and T. Combining this data for various Vin{= L — 2Nb + 1), Nb, we find that the relations 
(||)i(|l) describe the results quite well over a few orders in magnitude (Fig. ||). We have also included the data of Q| 
which used stochastic thermostats and see that the formulas work quite well. 
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FIG. 2. (left) AD/J^ against Vtn for Nb ^ 1 (□), Nb = 40 (o), with NH thermostats (left) and demons (right) at T = 0.5. 
Vi„/ (kX'^^^T^) (solid) and Vi„/(KA5^„j.r^)(l + 2a«:/Vi„) (dashes) are plotted. Application of the formulas to the data of [4] 
(A) works very well (data and formulas rescaled by 5000 for plotting). 




There are a few subtleties which we now resolve: First, we can and have made the distinction between the total 
volume and Vin in the formula (^, since we have performed analyses with Nb = 1 40 sites. Secondly, Cd can 
and does depend on T and also on the type of thermostats used, including Nb- Cd has similar behavior both for 
NH thermostats and for demons, with the former being larger, and both decreasing with Nb- Since the demons have 
twice the number of degrees of freedom per thermostatted site compared to NH thermostats, this is quite remarkable. 
Furthermore, somewhat surprisingly, when we increase Nb by one, thereby increasing the total number of degrees of 
freedom by 6(8) in the NH thermostat (demons) case, Cd decreases and so does AD for the same J. The reason for 
this will be clarified below. 

We study examples of the behavior of J2f=i ^^S- § (left). We see that the entropy production relation Eq. (^) 
works well near and far from equilibrium, and that the quadratic behavior with respect to J can be observed close 
to equilibrium. This is true, independently of the type of thermostats used. On the other hand, Af^^j^ can and does 
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depend on the type of thermostats used as well as T (Fig. || (right)). This is responsible for the thermostat dependence 
of C_D seen in Fig. ^. The larger A 
of more thermostats lead to more chaotic behavior 



J^jjj. for demons and the increase seen with Nb are quite natural since the existence 





and their quadratic 
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FIG. 3. (left) and j(l/rf - with respect to -J, for T = 0.5 (□,+), T = 2 (o, 

behavior Eq. (^ near equilibrium (solid) for the L = 162, Nb ~ 1 case. The quantities agree excellently. 

(right) The dependence of the X^ax on &t T = 0.5 for NH thermostats (□) and demons (o). This demonstrates that X'^^ 
is strongly heat-bath dependent, and is not uniquely determined by the local equilibrium conditions. 



We now study in detail the difference between AZ?, /^.D^in and lS.D„^ax'- Let us first work out what the condition 
AZ? < 1 means. Define r = AT/T, where T is the central temperature and AT is the difference in the boundary 
temperatures so that r < 2, always, and VT ~ rT /Vm- Then, we can derive a simple estimate, 
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AD<1 ^ v^^>^^^^r--^ (7) 

Let us analyze T = 0.5 case more concretely; the condition is most stringent for Nb — 1 (NH) case, when Xf^ax 
the smallest. So we obtain the condition Vin/r'^ ^ 130 for AD < 1 and for large lattices, Vin ^ 400, it should be 
satisfied for any gradient, which is consistent with our results. While we did not consider here the non-linearity of the 
profiles, the boundaries temperature jumps and the non-linearity of the response |jl^,|l8|, we are in some sense close 
to equilibrium when AD < 1 so that the rough arguments suffice for the purpose at hand. 

When AD < 1, how large is the difference ADmax — ^D? In this case, ADmin = AT), as explained above, so we need 
to consider \max + \nin , which is zero in equilibrium since the Lyapunov spectrum is symmetric with respect to sign 
inversion j^]. As we move away from equilibrium, the behavior can be described by Xmax/Xmin + 1 = CxJ^ + O {J^)- 
Then, when AZ? < 1, 

A^ "^^^-"^ 

For TVb = 1, T = 0.5, Cx = 13L°■^ so that {AD^ax - AD)/ AD < 0.1 x (y„/100)"° '5. Since the statistical errors 
are typically at a few % level, the difference between AD and AD^ax is at most comparable to them except for small 
systems, when AD < 1. In Fig. |l|, it can indeed be seen that for AD < 1, AD^ax agrees with AI? within error. 
Similar analysis can be applied at different T. 

When AD > 1, the relation ADmin — AD no longer holds exactly and the deviations from it is of 0{SX/X), 
with SX being the sum of the differences of the exponents Xj (j > Dky), which is not zero even in equilibrium. 
SX/X ^ AD/D so that in systems with small relative dimensional loss, ADmax can still be a good approximation. 
In Fig. 1^, we see that when AD > 1, AD^ax is in general not a good approximation to AD and that it is a clearly 
better approximation there for the larger system. 

We derived the extensivity of AD in systems with bulk behavior under thermal gradients in (^, thereby making 
the notion precise. The extensivity was explicitly related to the macroscopic transport properties as in d) through 
entropy production. Previously, it was emphasized that the extensivity of AD is not compatible with local equilibrium 
so that it is questionable |^ . It is now known that systems such as 0^ theory in = 1 — 3 display violations of local 
equilibrium under thermal gradients in a similar manner, as ~ (VT/T)^ [l8| . This resolves the apparent conflict since 
the violations of local equilibrium emerge in a manner analogous to (2). 

We further explicitly verified that AD in (f>^ theory behaves extensively under various thermal gradients for L ^ 10"*^ 
and its relation to transport. The relations (||),(§) and (||), however, are more general. We saw that the relations 
applied well to [Q which used stochastic thermostats in a different model. An application to dilute gas using the 
standard estimates of Xmax yields Co — 2/[pw^ ln(4^/d)], where p is the density, v is the average particle velocity, I 
is the mean free path and d is the particle diameter. Then for VT/T ~ 0.01 m~^, AD ^ 10* m~'^ at room temperature 
— quite large, yet far smaller than the total number of degrees of freedom. 
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We find the results satisfying from the physics point of view: Since AD pertains to the whole system, it includes the 
temperature profile which is curved in general, boundary temperature jumps and the various types of thermostats. 
Yet, AZ? can be related to macroscopic transport with the thermostat dependence cleanly separated into Xmax- 
Furthermore, AD has extensive behavior with respect to the internal volume wherein the system is manifestly in 
non-equilibrium. We have seen that A^^^. is not unique: In global thermal equilibrium, different choices of heat-baths 
can lead to very different values. The result is that dimensional loss is not unique either, only the product ADX'^^^ 
behaves macroscopically and can be related to thermodynamic quantities. 
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